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Abstract 

We propose two new definitions of tlie exponential function on 
time scales. The first definition is based on the Cayley transformation 
while the second one is a natural extension of exact discretizations. 
Our eponential functions map the imaginary axis into the unit cir- 
cle. Therefore, it is possible to define hyperbolic and trigonometric 
functions on time scales in a standard way. The resulting functions 
preserve most of the qualitative properties of the corresponding con- 
tinuous functions. In particular, Pythagorean trigonometric identities 
hold exactly on any time scale. Dynamic equations satisfied by Cayley- 
motivated functions have a natural similarity to the corresponding 
diferential equations. The exact discretization is less convenient as far 
as dynamic equations and differentiation is concerned. 
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1 Introduction 



The main goal of the time scales approach is a unification of the differential 
and difference calculus [H |2l [3]. There are many papers and books where 
notions and notation concerning time scales are explained in detail [H El [5], 
see also [1]. A time scale T is defined as an arbitrary closed subset of M [Il[3]. 
The forward jump operator a is defined as a{t) = inf {s G T : s > t} (and we 
assume inf = sup T). We usually denote t" := crit). 

A point t G T is called right-dense iff t*^ = t and right-scattered iff 
> t. Similarly one can define the backward jump-operator p, left-dense 
points and left-scattered points [T]. Sometimes we have to exclude from 
considerations the left-scattered maximum of T (if it exists). T minus the 
left-scattered maximum will be denoted by T*^ (if such maximum does not 
exist then T'^ = T). 

Graininess p = fi{t) is defined, for t G T'', as /i(t) = t" — t. The delta 
derivative (an analogue of the derivative for functions defined on time scales) 
is defined by 

.Ht):- lim ^K;))-^M . ((,t«). (1.1) 

S t CF[t) - S 

If t is right-dense, then x^{t) = x{t). If t is right-scattered, then is the 
corresponding difference quotient. 

An important notion is rd-continuity. A function is said to be righ-dense 
continuous (rd-continuous), if it is continuous at right-dense points and at 
left-dense points there exist finite left-hand limits. The graininess /x is rd- 
continuous but, in general, is not continuous. Dynamic equations are time 
scales analogues of differential equations (i.e., they may contain delta deriva- 
tives, jump operators and sometimes the graininess of the considered time 
scale) [U [5] . 

In this paper we consider the problem of defining elementary (or special) 
functions on time scales. We propose a new approach to exponential, hyper- 
bolic and trigonometric functions. We point out that the definition of such 
functions is not unique. We suggest new definitions (improved exponential, 
hyperbolic and trigonometric functions) based on the Cayley transformation. 
The new functions preserve more properties of their continuous counterparts 
in comparison to the existing definitions. In particular, our exponential func- 
tion maps the imaginary axis into the unit circle and trigonometric functions 
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satisfy the Pythagorean identity. Dynamic equations satisfied by our im- 
proved functions have a natural similarity to the corresponding differential 
equations. 

We also propose the notion of exact time scales analogues of special func- 
tions and we identify Hilger's definition of trigonometric functions |6] with 
exact analogues of these functions. Exact discretizations of differential equa- 
tions [3 IHl E] are intended as a way to connect smooth and discrete case, but 
in an apparently different way than the time scales calculus. In this paper 
we incorporate the exact discretizations in the framework of the time scales 
approach. We discuss exact exponential, hyperbolic and trigonometric func- 
tions on constant time scales. An extension of these results on more general 
time scales seems to be rather difficult. 

2 Survey of existing definitions 

The exponential function on time scales has been introduced by Hilger [T], 
and his definition seems to be commonly accepted [21 IS] • A different situation 
has place in the case of hyperbolic and trigonometric functions, where Hilger's 
approach [6] differs from Bohner-Peterson's approach pU]. 

2.1 Exponential function 

Hilger defined the exponential function as follows: 

ft 



ea{t, t) ■= exp (^J ^^(,)(a(s))As^ 



e.W :=e„(t,0) , (2.1) 



where 



a(^) := y log(l + zh) (for h>0) and Co{z) := z . (2.2) 
n 

This definition applies to the so called ^-regressive functions a = a{t), i.e., 
those satisfying 

1 + ^(t)a(t) ^ for all teJ'" . (2.3) 

Such functions are usually called regressive, but we reserve this name for 
another class of functions, see Definition 13.11 
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In the constant discrete case (T = eZ, a = const) we have 



e„(t) = + , (2.4) 
and in the case T = R we have 

ea{t) = exp / a{r)dT . (2.5) 

^0 

Theorem 2.1 ([Ulin])- // a, /9 : T — )■ C are ^-regressive and rd- continuous, 
then the following properties hold: 

1. ea(r,to) = (1 + A^(t)a(t)) e„(t,to) , 

2. {ea{t,to)) ^ = eeMa(t,to) , 

3. eait,to) eaito,ti) = ea{t,ti) , 

4- Gait,to) ei3{t,tQ) = Ca^i^isitjto) , 

where a ©'^ /3 := a + /? + and Q'^a 



The addition ©'^ is usually denoted by ©. However, we reserve the no- 
tation © for another addition, see Definition 13.81 The exponential function 
(12. ip solves the Cauchy problem: = ax, x(0) = 1, see p. 

Similar considerations (with the delta derivative replaced by the nabla 
derivative) lead to the nabla exponential function |Tl]. In the case T = 
the nabla exponential function is given by 

e„(t) = (1 - ae)-^ (2.6) 

and solves the Cauchy problem: = ax, x{0) = 1. The existence of the 
above two definitions refiects the duality between delta and nabla calculus. A 
linear combination of delta nad nabla derivatives (the diamond-a derivative 
[12] ) leads to another definition of the exponential function (the so-called 
diamond-alpha exponential function [T3]). 
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2.2 Hilger's approach to hyperbolic and trigonometric 
functions 

The first approach to hyperbohc functions has been proposed by Hilger [B], 

cosh.(t) = '-^'^ , sinh„(t) = ^"(^^ -^'^'-^'^ , (2.7) 

where a is yU-regressive. Among its advantages we have the identity 

cosh^(t) -sinh2(t) = 1 . (2.8) 

Delta derivatives of these functions are linear combinations of both hyperbolic 
functions, e.g., 

cosh^(t) = ^ ^ cosh,(t) + ^ ^ sinh„(t) . (2.9) 

In the constant discrete case (T = eZ, a = const) we have 
'1 + ae)e + (1 + ae)-e 



cosho,(t) 



sinhQ,(t) 



2 

''1 + ae)^ — (1 + as) 



(2.10) 



The hyperbolic functions (12. 7p evaluated at imaginary a are not real. Thus 
the definition (12. 7p cannot be extended on trigonometric functions by substi- 
tuting iu for a. In order to treat trigonometric functions Hilger introduced 
another map a — j- lu, see [6]. In the constant case (i.e., n{t) = const and 
u = const) the final result is very simple 

coSaj(t) = coscjt , sin(^(t) = smut . (2-11) 

In fact, this is a particular case of the exact discretization (see Section H]), 
although Hilger's motivation seems to be different. Exact discretizations have 
many advantages (compare [14J) but, unfortunatelly, their delta derivatives 
are quite complicated. For instance, at right scattered points we have: 

A sinoj/i cosw/i — 1 

(smojt) = cosut-l smut, (2.12) 

/i n 

compare [B], see also Section HI 
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2.3 Bohner-Peterson's approach to hyperbolic and tri- 
gonometric functions 

The second approach has been proposed by Bohner and Peterson (TUl IS] 
cosh„(t) = , 

(2.13) 

sinh„(f) = , 

where a is /i-regressive. The hyperbohc functions defined by fl2.13p satisfy 

cosh^(t) = asmh.a{t) , sinh^(t) = acoshQ,(t) , (2-14) 

The identity (12. 8p is not vahd. Instead we have 

cosh^(t) - sinh2(t) = e„^,2(t) . (2.15) 

Bohner and Peterson define trigonometric functions in a natural way, evalu- 
ating hyperbolic functions at the imaginary axis 

cos^(t) = coshi^(t) , isin^(t) = sinhi^(t) . (2.16) 

Then, 

cos^(t) + sin^(t) = e^^2(t) . (2.17) 
In the constant discrete case (T = eZ, a = const, u = const) we have 



(1 + «£:)=+ (1 - ae)^ 
coshQ(t) = 



sinhQ,(t) 
Moreover, 



2 

[1 + ae)^ — (1 — tt£) 



(2.18) 



e_,,2(t) = (1 - e^a^ , e,^2(t) = (1 + eV)- . (2.19) 
Therefore, 

lim e_£Q2(t) = oo , lim e^sa^it) = , (2.20) 

t— >— oo t—^oo 

provided that \ae\ < 1. Similarly, 

lim e^^2{t) = , lim e£^2(t) = oo , (2.21) 

t—^—oo t— >oo 

Therefore, the definition fl2.13p leads to Pythagorean- like identities fl2.15p . 
f l2.17p which have essentially different behaviour in the discrete and contin- 
uous case. 
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3 Approach motivated by the Cayley trans- 
formation 



In this section we present new definitions of exponential, hyperbolic and 
trigonometric functions and their properties. We will tentatively refer to 
them as 'improved' functions, because they simulate the behaviour of their 
continuous counterparts better than the previous definitions. Our new defi- 
nitions are based on the classical Cayley transformation: 

z ^ caj{z,a) = Y^-^ , (3.1) 

see, for instance, [16] . 



3.1 New definition of the exponential function 

In order to formulate our definition we need to redefine a notion of regres- 
sivity. 

Definition 3.1. The function a : T E C is regressive if fj,{t)a{t) ^ ±2 for 
any t eT'^. 

Definition 3.2. The improved exponential function (or the Cayley- exponen- 
tial function) on a time scale is defined by 

E^{t,to) := exp (^^ Cm(s)(«(s))As^ , E^{t) := E^{t,0) , (3.2) 

where a = a{t) is a given rd- continuous regressive function and 

a(^):=-^log[^|^ (for /i>0) andCo(^):=^. (3.3) 

Here and in what follows the logarithm is understood as a principal branch 
of the complex logaritm with image [—iir, in]. 

Lemma 3.3. If a is rd-continuous and regressive, then the delta-integral in 
( f5^) exists. 
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Proof: The assumption of regressivity of a implies that the logarithm in (j3.2p exists 
(is finite) for any t € T*^. Thus the function t — >■ ^^(t)(a(t)) has no singularities. 
To complete the proof we will show that o a is rd-continuous (which implies 
that it has an antiderivative, see [1]). At right-dense to we have 

lim a{t) = a{to) , lim n{t) = /u(to) = , (3-4) 
because a and fi are continuous at right-dense points. Therefore 

lim C^,{t){a{t)) = hm — log f = a{to) . 3.5) 

i-i-io ' t^U) fi{t) 1 - ^a[t)fi{t) 

On the other hand, C/i(t(,)(a(io)) = Co(o(io)) = a(io)- Therefore, C^tO" is continuous 
at right-dense points. At left-dense sq we denote 

a(sQ ) := lim a{t) , ^(^0 ) ■= 1™ ^(0 = ■ (3-6) 

In general a(so ) 7^ ^(so), m(so ) 7^ ^(■so) but rd-continuity guarantees that all 
these values are finite. Then, 

hm C^,(t){a{t)) = hm -ttt log = a{sQ ) , (3.7) 

and the existence of this finite limit means that QnO a is rd-continuous. □ 



In the constant discrete case (T = eZ, a = const) we have 

Eait) = (\^) ^ , (3.8) 

and in the case T = R we have, as usual, 

Ea{t) = exp [ a{T)dT . (3.9) 
Jo 

The formula (13.81) appeared earlier in different contexts, see for instance 
[Ml [29]. 

The new definition (13. 2 p of the exponential function can be related to 
Hilger's definition (12. ip with an other exponent. Namely, we are going to 
prove that Eait^to) = epit^to) provided that 

m = , r% • (3.10) 



8 



Theorem 3.4. For any regressive, rd-continuous a = a{t) there corresponds 
a unique ^-regressive, rd-continuous 13 = Pit) (given by such that 

Ea(t,to) = epit^to). For fi-regressive, rd-continuous j3 satisfying ^ —2 
there exists a unique regressive, rd-continuous a given by 

= 1 ' (3-11) 

such that Ea(t,to) = ep{t,to). 



Proof: Ea{t, Iq) = ep{t^ to) if and only if the integrands in ()2.ip and ()3.2p coincide, 
i.e., i^{l3) = C^t(a)- Thus /3(t) = a{t) at right-dense points, and for / (i.e., at 
right-scattered points) we get 

1 -I- iita 

l + /i/3 = - f^. (3.12) 

1 — 2^*0 

Both cases lead to a single condition (jS.lOp (or, equivalently, to (jS.lip ). To com- 
plete the proof we verify that ///3 = —1 iff = —2. Therefore for any regressive 
a we have ///3 / —1. Then, = 2 corresponds to /i/3 = ±oo. Thus for any 
^-regressive (3 we have |/ia| ^ 2. Finally, we observe that /i/3 = —2 corresponds to 
fia = zbcxD, and for any other values of /3 the value of a is uniquely determined by 

dsn]). 

To complete the proof we have to show that a is rd-continuous iff /3 is rd- 
continuous. At right-dense to functions a, fj, are continuous and also /i(to) = 0. 
Therefore, from (13. lip we get 

lim a{t) = lim /3(t) , a{to) = (3{to) . (3.13) 

Hence, a is continuous at to if and only if /3 is continuous at to- If so is left-dense, 
then //(t) — 7> as t — )• Sq . As a consequence, we have 

a(so ) = lim_ a{t) = lim_ /3(t) = /3(so ) . (3.14) 
Therefore, ol{sq) exists if and only if /3(sq ) exists, which ends the proof. □ 



3.2 Properties of the Cay ley- exponential function 

First of all, we observe a close relation between C,^ and the Cayley transfor- 
mation (13. ip . 

e'^fft(^) = cay(z, ^/i) . (3.15) 
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The transformation inverse to is given by 

^ = C(C) = ^tanh^ (MO), Co-'(C) = C. (3.16) 
In particular, 

2;=— tan^ for Q = iu . (3-17) 

it ^ 

Therefore, maps the open segment (— X'x) ^ ^■'^ onto zR, and M is 
mapped onto the real segment (~f > f )■ 

Corollary 3.5. C/i maps the imaginary axis onto the segment (— ^, ^) C iM. 



Lemma 3.6. Denoting ( = 'y + ir] and taking into account Ii3.16\) . we have 
l^l < — <^=^ cos?7/i > , 



h 
2 

h 



cos Tjh = , (3.18) 



, , 2 

\z\ > — <^=^ cosr]h < . 
h 

Proof: We compute: 

T ~ e'J7+*/*r? + 1 ~ e^'y COS /it/ + 1 + ie^T sin hrj ' ^ ' 



tanh — 
2 



' e^^^T - 2e'*T cos /it/ + 1 



(3.20) 



y e^'^T + 2e''i' cos hr] + I ' 
To complete the proof it is enough to notice that e'*'*' > and ^ = tanh ^ . □ 

Corollary 3.7. (h maps the disc \z\ < onto the strip < ^ < ^■ 



10 



Definition 3.8. Given a time scale T and two functions a, /3 : T — >■ C, we 

define 

where ji = ii{t) is the graininess ofT. 

Lemma 3.9. The function has the following properties: 



Cm(«) = ' C/u(-Q;) = -Cm(«)> C-/u(a) = Cm(«) > (3-22) 

C(«) + Cm(/5) = Cm(«©/3), (3.23) 
where bar denotes the complex conjugate and we assume (in order to avoid 
infinities) ^/la ^ —1, ^ —1 and ^ — (|//Q;) ^ 

Proof: The function C,^ = Cuict) is analytic with respect to a (and n is real). Hence 
the first property follows. Other properties can be shown by direct calculation: 

1 1 - Ian 1 1 + ^afi 

Cm -« = - log T— T = -- log 1 = -C^. a , 

M 1 + ^a/j, n 1 — ^a/i 

C-m(«) = -- log — f-^ = - log = Cm(«) , (3.24) 

/i 1 + ^a/x /i 1 — ^cKyu 



Ua) + CM = - log .f,. r = CM © 

1 + j/x^ap - |/i(Q; + p) 



provided that jfi'^aP / —1. □ 



Theorem 3.10. // a, /3 : T — >■ C are regressive and rd- continuous, then the 
following properties hold: 



1 + ^ljL{t)a{t) 
2. (EMto))-^ = E_Mto) , 



3. Ea{t,to) — Ea{t,to) , 
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5. Ea{t,to) Ei3{t,to) = Ea^l3it,to) , 

where we use a standard notation t" = (j{t). 

Proof: It is sufficient to prove the first property for right-scattered points {f^ > t). 
E^if'M) = exp ^ ^ C^(,)(a(r))Arj exp (^^ C^(,)(a(r))AT^ (3.25) 



Then, using (j3.3p . we get 

1 + -u(t)a(t) 

f C,irMr))Ar = /.(t)CMt)(«(i)) = log ^ f)' \{ , 
t 1 - ^fi(t)a{t) 

and substituting it into (I3.25P we get the first property. The second property 
follows directly from C,^{—a) = —Cfi{oi), see Lemma Indeed, 

Ea^it,to) = exp i- C^i(T){a{r))AT] = exp CA»(T)(-a(r))Ar = to). 
V Jto / Jto 

The third property follows directly from analycity of the exponential function and 
from Lemma 13.91 We recall that t € T C M. Indeed, 

Ea(.t,to) = exp / Cf.{r){a{r))AT = exp / Cf^(.T)(M^)AT = E^{t,to) . 
Jto Jto 

The fourth property is derived in a straightforward way: 

Ea{t,to) Ec,{to,ti) = exp ( [ CMt)(«W)At+ [ \f,(^t){a{t))At) = E^{tM)- 

\Jto Jti J 

Finally, 



E^(tM) Efi(tM) = exp / (Cm*) («(*))+ Cm*) At = S„e^(t,to), 

Jto 

where we took into account Lemma |3. 91 □ 



The formula (13.211) is identical with the Lorentz velocity transformation 
of the special theory of relativity (the role of the speed of light c is played 
by ^). Denoting 



a := -/io; , j3' = -///? . (3.26) 
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we can rewrite the formula fl3.21l) in a simpler form 

"'®^'-T^. (3.27) 



Lemma 3.11. If a' and (5' are real functions on T and a' © is given by 
(d^, then 



\a'\<l and |^'| < 1 =^ |a'©/3'|<l, (3.28) 

\a' ©^'1 = 1 \a'\ = 1 or \f3'\ = 1 . (3.29) 

Proof: Using (j3.27p we compute 

which immediately yields (j3.29p . Then, \a'\ < 1 and |/3'| < 1 imply that the righ- 
hand side is positive. Hence, we have (|3.28p . □ 



Definition 3.12. The function a : T — )■ M zs called positively regressive if 
for all t G T'' we have |a(t)/i(t)| < 2. 

Theorem 3.13. // a : T — > M rd-continuous and positively regressive, 
then the exponential function is positive (i.e., Ea{t) > for all t G T). 



Proof: If a is real and positively regressive, then (for any t G T*^) we have 

> . (3.31) 



1 + ^/i(t)a(t) 



1 - \^i{t)a{t) 

Thus C^(j)(a(t)) is real for any t G T'' and, as a consequence, the exponential 
function is positive. □ 

The regressivity condition \a{t)^{t)\ < 2 is automatically satisfied at 
right-dense points. At right-scattered points we have 

E.if^) - E^{t) = ^{t)a{t)^^^^^^^^ . (3.32) 

13 



Therefore, the condition \a(t)ii{t)\ < 2 is equivalent to 

\E^{r) - E"{t)\ < \E^{n + E^{t)\ (3.33) 
In the real case it means that Ea{t'^) and Ea{t) have the same sign. 

Theorem 3.14. The set of positively regressive real functions is an 
abelian group with respect to the addition ©. 

Proof: The formula (j3.2ip obviously yields a © /3 = /3 © a. The element inverse to 
a (given simply by ©a = —a) always exists in TZ^. Therefore, it is sufficient to 
show that TZ^ is closed with respect to ©. Taking into account (I3.2ip we have: 

l + ^fia(Bf3 ^ (1 + ^^xa){l + 
l-^fxa®f3 (1 - i/ia)(l - 

(compare the last line of (|3.24p ). If a,/3 G 7^+, then all four terms at the right- 
hand side are positive. Therefore, the left-hand side is positive which means that 
\^ia®f3\ < 2. □ 

The set 71 of all regressive functions is not closed with respect to ©. In- 
deed, suppose that a regressive a is given. The formula fj,'^a(3 = —4 uniquely 
determines (3 which also has to be regressive. However, in this case a (B f3 
becomes infinite. 

Lemma 3.15. 

x^{t) = I3{t)x{t) ^ x^{t) = a{t){x{t)) , (3.35) 
where j3 is given by Ii3.10\) and 

(xit)) := ^^'^ . (3.36) 

Proof: Direct computation shows 

x^ = a{x) ^ 2x^{l + ^fiP) = Pix + x") ^ x^ = /3x, 
where first we used p. lip , and then we substituted x'^ = x + jix^ . □ 
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Theorem 3.16. The exponential function Ea{t,tQ), defined by liS. 2\] . is a 
unique solution of the following Cauchy problem: 

x^{t) = a{t){x{t)) , x{to) = 1 , (3.37) 

where a is regressive rd-continuous function and {x{t)) is defined 1^3. 36\) . 



Proof: By Lemma 13.151 the initial value problem (j3.37p is equivalent to 

x^{t) = l3{t)x{t) , x{to) = 1 . (3.38) 

Taking into account that /i and a are rd-continuous, we have /? rd-continuous 
(compare Theorem 13. 4p . Therefore, we can use Hilger's theorem concerning the 
problem (|3.38|) . stating that its unique solution is efs{t,to). Finally, we use Theo- 
rem [331 once more. □ 



In the discrete case the equation fl3.37p (treated as a numerical scheme) 
can be interpreted either as the trapezoidal rule, implicit midpoint rule or 
the discrete gradient method [T71 HI |T9] . These implicit methods are more 
accurate than the explicit (forward) Euler scheme (which is related to the 
equation x^ = ax) and can preserve more qualitative, geometrical and phys- 
ical characteristics of the considered differential equations, e.g., integrals of 
motion ^19j . 

3.3 New definitions of hyperbolic and trigonometric 
functions 

Our definition of improved hyperbolic and trigonometric functions follows in 
a natural way (similarly as in the continuous case) from the definition of the 
exponential function. 

Definition 3.17. Cayley-hyperbolic functions on a time scale are defined by 
Cosh„W := M)±^ , Si„h„W := ^^W^^ , (3.39) 

where the exponential function is defined by jiS. 2\] . 
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Definition 3.18. Cayley-trigonometric functions on a time scale are defined 
by 

Cos„W := ^^^^^^ . Sin„(i) := . (3.40) 

In other words, Cos^it) = Cosh.i^{t), zSin^(t) = Sinhj^(t). 

Properties of our hyperbolic and trigonometric functions are identical, or 
almost identical, as in the continuous case. Note that below we often use the 
notation fl3.36p . 

Theorem 3.19. Cayley-hyperbolic functions satisfy 

Cosh^ (t) - Sinh^ (t) = 1 , (3.41) 



Coshf (t) = a(t)(Sinh„(t)) , Sinhf (t) = a(t)(Cosh„(t)) . (3.42) 

Proof: By Theorem 13.101 we have Ea{t)E-a{t) = 1. This is sufficient to directly 
verify the identity (|3.4ip . By Theorem 13.161 we have E^{t) = a{t){Ea{t)) and 
differentiating IK39\i we get ([332]). □ 



Theorem 3.20. Cayley-trigonometric functions are real-valued for real u 
and satisfy 

Cos^(t) + Sin^ (t) = 1 , (3.43) 



Cos^(t) = -u{t){Sm^{t)) , Sin^(t) = a;(t)(Cos^(t)) . (3.44) 



Proof: By Theorem I3.1UI we have Eicj{t) = E^i^{t). Thus the reahty of trigono- 
metric functions fohows directly from Definition 13.181 Using Ei^{t)E-i^{t) = 1 we 
get the Pythagorean identity (|3.43p starting from ()3.40p . Derivatives ()3.44p can 
be obtained by straightforward differentiation of the exponential functions using 
Theorem 13.161 □ 

Theorem 3.21. The function C 3 a Ea{t) maps the imaginary axis into 
the unit circle, i.e., Iiea{t) = =^ \Ea{t)\ = 1. 
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Proof: We compute 

\Ea{t)\^ = E^{t)Ejt) = Ea{t)E^{t) = E^^a , (3.45) 

where we used twice Theorem I3.10[ From the formula ()3.2ip we see immediately 
that 

a © a = -4=^ a = —a , 
i.e., a = ioj, where uj{t) € M for t G T. Therefore, Ea^ait) = 1 for imaginary a. □ 

Theorem 13.211 is crucial for a proper definition of trigonometric functions. 
Indeed, \Ei^\ = 1 is equivalent to (ReEi^Y + (ImEi^Y = 1 and we can 
identify HeEi^^ and IraEi^^ with trigonometric functions. 

At the end of this section we present second order dynamic equations 
satisfied by hyperbolic and trigonometric functions with constant a and cu, 
respectively. 

Lemma 3.22. Averaging commutes with delta differentiation, i.e., 

{x{t))^ = {xHt)) . (3.46) 

Proof (x(t))^ = i {x{t) + xit^f = \ {x^{t) + x^in) = (x^(t)) . □ 

Proposition 3.23. // a{t) = const, then improved hyperbolic functions on 
a time scale, Cosho, and Sinha, satisfy the equation 

x^^(t) = a^{x{t))) , (3.47) 

where {{x{t))) = \{x + x" + x'''') and x"" {t) := xit") . 

Proof: follows immediately from (j3.42p and (|3.46p . Indeed, 

Cosh^^(t) = a(Sinh„(t))^ = a(a(Cosh„(t))) = a2((Cosh^(t))) . 
The same calculation can be done for Sinh^^. □ 

Proposition 3.24. // u{t) = const, then improved trigonometric functions 
on a time scale, CoS(^ and Sin^, satisfy the equation 

x^^ + a;2((x(t))) =0 . (3.48) 
Proof: by straightforward calculation, compare the proof of Proposition 13.231 □ 
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4 Exact special functions on time scales 



The simplest (almost trivial) way to construct time scales analogues of special 
functions is to take their exact values. 

Definition 4.1. Given a function f : M. ^ C we define its exact analogue 



Although the path / — )■ / is obvious and unique, the inverse correspon- 
dence may cause serious problems. Usually, it is not easy to find or to indicate 
the most appropriate (or most natural) real function corresponding to a given 
function on T. 

4.1 Exact exponential function on time scales 

The continuous exponential function is given by ea{t) = exp J^* a{T)dT, where 
a : R — )■ C is a given function. A non-trivial question is to choose a function 
a = a{t), provided that we intend to define an exact exponential correspond- 
ing to a given function a : T — > C. In general, the choice seems to be highly 
non-unique. 

In this paper wc confine ourselves to the simplest case, a — const, when 
it is natural to take a{t) — a — const. 

Definition 4.2. The exact exponential function E^{t, to) (where a is a com- 
plex constant and t,to eT) is defined by E^{t,to) :— e'^(*~*o). 

Theorem 4.3. The exact exponential function E^^{t,to) satisfies the dy- 
namic equation 



f :T as f := /|t, i.e., 



fit) := fit) (for teT). 



(4.1) 



aipait){xit)) , 



xito) = 1 , 



(4.2) 



where ipait) 



1 for right- dense points and 



anit) 



2 



tanh 



a^it) 
2 



(4.3) 



for right- scattered points. 
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Proof: For right-dense t the equation (j4.2p reduces to = ax. Then, still 
assuming t right-dense, we compute 

^ ' dt 

i.e., e"^*"*"-* satisfies the equation x^ = ax. For right-scattered t we have: 
E^J^if'^to) - E^^{t,tQ) = e"(*"-*o) - e"(*-*o) = e"(*-*o) (e"'^ - 1) , 
E^'^if^to) + E^^it^to) = e°(*"-*o) + e"(*-*o) = e"(*-*o) (e"'^ + 1) , 

and, therefore, 

S-(r,to) - i?r(i,io) = tanh^ (i?-(r,to) + i?r(i,io)) , (4.4) 



which is equivalent to (j4.2p . The initial condition is obviously satisfied. □ 

The equation (14. 2 p can be interpreted as the dynamic equation fl3.37p 
with a modified a (i.e., a — )■ ailja). Another interpretation can be obtained 
by a modification of the delta derivative. We define 

:,A'„(^) lim "^f,^ ~ "^^"^ (4.5) 

where 5a is a function given by 

5q(/x) := — tanh-^ . (4.6) 
a 2 

Lemma 4.4. 

x^(t) = ^«(f)x^Hi) (4.7) 



Proof: We compute 

X(t'^)-X(s) X(t'^)-X(s) dJf-s) A' , X , ,^ 

s^a(t) s^a(t) s=ia{t) 

which yields (gZ])- D 
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Corollary 4.5. The equation (^T^, satisfied by the exact exponential func- 
tion, can be rewritten as 

= a{x{t)) , x{to) = 1 . (4.8) 

The equation (14.21) is the exact discretization of the equation 
In general, by the exact discretization of an ordinary differential equation 
X = f{x), where x{t) G M^, we mean the difference equation Xn+i = F{Xn), 
where X„ G M^, such that X„ = x(t„). Any equation has an implicit exact 
discretization (provided that the solution exists). It is worthwhile to point 
out that all linear ordinary differential equations with constant coefficients 
have explicit exact discretizations [3 El [9]. 



4.2 Exact hyperbolic and trigonometric functions on 
time scales 

In order to simplify notation we confine ourselves to to = 0. All results can 
be obviously extended on the general case. 

Definition 4.6. Given a real constant a, exact hyperbolic functions on a 
time scale T are defined by 

cosh^'^(t) = — — — — = cosh at , 

(4.9) 

sinhr(t) = ^^^^^^^^^ = sinhat. 
2 

We point out that in this definition t G T. The same remark concerns 
the next definition. 

Definition 4.7. Given a real constant u, exact trigonometric functions on 
a time scale T are defined by 

, EZ{t) + EtUt) 

COS^ it) = = COS Ujt , 

^ (4.10) 

. , Ef^^jt) - EZUt) . . 
sm^ (t) = — = sm cut . 

It turns out that dynamic equations satisfied by exact trigonometric and 
hyperbolic functions are rather awkward in the case of arbitrary time scales. 
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These equations simplify considerably under asumption of constant graini- 
ness. Therefore, from now on, we confine ourselves to the case = const. 

Moreover, we observe that the function -0^ (defined by (14.31) ) is symmetric. 
In particular, i/jiuj(t) = i/j-i^it). We denote: 

(f){x) := - tan ^ (for x ^ 0) , (/.(O) = 1 . (4.11) 
Therefore, for fi{t) = const we have 

ipiuiif) = (pi^fJ-) = const . (4-12) 

Theorem 4.8. We assume that the time scale T has constant graininess ^. 
Then exact trigonometric functions on T, i.e., cos^^ and sin^^, satisfy the 
dynamic equation 

x^^(t) + uj^(t)\uj^i) {{x{t))) = , (4.13) 
which is equivalent to 

c^'^it) + (sinc^^) ' x(r) = , (4.14) 



X 



where 



sin X 

sinc(x) := (for x ^ 0) , sinc(O) := 1 . (4-15) 



X 



Proof: By Theorem 13. 161 we have 

{E-^,J^ = ±iu:^^Ut){Ef^J . (4.16) 



Therefore, taking also (j4.1ip into account, we get 

(cos-(t))^ = i {Et^{t) + EtUt)f = -a;0(a;M)(sin-(t)) 
(sin-(i))^ = i {EZit) - E-_Ut)f = ^<A(a;/x)(cos-(t)) . 



(4.17) 



We notice that the second formula is equivalent to (j2.12p . Applying Lemma 13.221 
to (|4.17|) . we obtain 

(cos-(t))^^ = -u^^<p\u;f,){{costrm , 

(4.18) 

(sin-(t))^'^ = -u;^^\uf^){{smtfm , 
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i.e., we have (j4.13p . The equivalence between ()4.13p and ()4.14p is obvious for 
/i = 0. In the case /x 7^ we use trigonometric identities: 

sin(cjt) + sin(a;t + 2a;^) + 2sm{ujt + uj/i) = 2(1 + cosw/u) sin(a;t + w/u), 

(4.19) 

cos(a;t) + cos(a;t + luj/i) + 2cos(Li;t + ufi) = 2(1 + cosw/u) cos(a;t + ujfi), 
obtaining 

((sinS'(t))) = (cos^) sm{ujt + ujn) , 

{{cosff{t))) = (^cos'^^ cos{ujt + ujfi) . 

Taking into account that sinc^ = 4){ojh) cos ^ and = t + we get 

(l^-'iufiXixm = (sinc'^yxin , (4.21) 

where x{t) is an arbitrary hnear combination of sin^^ and cos^^, which ends the 
proof. □ 

The exact discretization of the harmonic oscillator equation leads to an- 
other modification of the delta derivative (see [T^ 122]). 

x^-(t) = lim ^ ^ ^ ^ ; ^ ^ ^ . 4.22 

s^t OJ-^ sm^uf" - us) 

In order to avoid infinite values of we assume \u^{t)\ < n. All posi- 
tively regressive constant functions u (see Definition I3.12p obviously satisfy 
this requirement. 

Proposition 4.9. If x = x{t) satisfies x + u'^x = for t eM., then 

(x(t)|ieT)^-=i(t)keT. (4.23) 

Proof: By assumption, x(t) = Acosut + B smut. Then 

x{t'^) = cos ufj, (A cos ut + B sin ut) + sin u^i {B cos ut — A sin ut) , 

because = t + fi. By direct computation we verify 

x{t") — x{s) cos{ut'^ — us) = {B cosws — ^sinws) sin(Li;t'^ — us) . 
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Therefore, 

A//,, ,. xif^) — x(s) cosiuf^ — UJS) ,^ , . , .,, 

x^-'(t)= hm ^ ^ , i ^ ^ -= uj(B cos ujt - Asmujt) = x(t), 

which ends the proof. □ 



Lemma 4.10. 



x^{t) = smc{ufi) - ^fiuj^ (^sinc'^y x{t) . (4.24) 



Proof: Substituting the definition (|4.22p into the right-hand side of the formula 
(I4.24p . and taking into account that 

sm(ujt^ — ujs) 
lim = smclujfi) , 

S tOf" - UJS 

hm = -fiuj I smc— 1 x{t) , 

s— >-i I — s z V z y 

we obtain the left-hand side of ()4.24p . compare (jl.ip . □ 



Proposition 4.11. The equation ^.l^^ , satisfied by exact trigonometric 
functions, can be rewritten as 



X 



A" A" 



{t) + u^x{t) = . 



(4.26) 



Proof: The formula (j4.24p can be rewritten as 

x'^ (t) = sinewy (^cos x^" (t) - w sin ^ x(t)) , (4.27) 
and, therefore, 

x^^it) = sinc^"^ (cos^ x^'^'^Ht) - u^isinuf,) x^^(t) + oj^ sin^ ^ x(t)) . 
Moreover, x{t'^) = x{t) + fix^{t). Hence 

x{r) = (cosw^) x{t) + x'^''^{t) . (4.28) 
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Thus we can easily verify that 



X it) +a; smc — x(t) = sine — 




±{^x^'^K(^t)+u^x{t)), (4.29) 



which ends the proof. 



□ 



The exact discretization of the harmonic oscillator equation x + u'^x = 
was discussed in detail in [H]. In particular, we presented there the discrete 
version of the formula fl4.26p and related results. 

5 Conclusions and future directions 

We proposed two different new approaches to the construction of exponential, 
hyperbolic and trigonometric functions on time scales. 

The resulting functions preserve most of the qualitative properties of the 
corresponding continuous functions. In particular, Pythagorean trigonomet- 
ric identities hold exactly on any time scale. Dynamic equations satisfied by 
Cayley-motivated functions have a natural similarity to the corresponding 
diferential equations. 

The first approach is based on the Cayley transformation. It has impor- 
tant advantages because simulates better the behaviour of the exponential 
function, both qualitatively (e.g., the Cayley-exponential function maps the 
imaginary axis into the unit circle) and quantitatively, because 



i.e., this factor approximates exp(a/i) up to second order terms, while l + fia 
or (1 — fia)~^ are only first order approximations. 

Our approach has some disadvantages, as well. The promising notion of 
complete delta differentiability [201 EI] becomes very difficult or impossible to 
apply (because integral curves of our new dynamic systems on time scales, like 
fl3.37p . are not completely delta differentiable). Moreover, dynamic equation 
on time scales become implicit (and equations of standard delta calculus on 
time scales are explicit). However, nabla calculus is also implicit, and it is 
very well known that implicit numerical finite difference schemes have better 
properties than explicit schemes, see for instance fL9\ . 

The second approach consists in exact discretization. In this paper we 
confined ourselves to the simplest case, i.e., to the exponential function Ea{t) 



1 + 



1 + a/i + -{afif + . . . 



(5.1) 



1 — 
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with a = const. It would be interesting to define and study exact exponential 
functions for a larger class of functions a. We leave it as an interesting 
open problem. Other problems are associated with finding dynamic systems 
which correspond to exact discretizations. In the case of linear equations this 
subject is well known, but more general results are diffcult to be obtained. 

We point out that definitions of elementary functions on time scales are 
not unique. We presented some arguments in advantage of our definitions 
but in principle one can develop several different theories of elementary and 
special functions on time scales. It seems important to develop and under- 
stand different approaches, and, if possible, to find a "vocabulary" to trans- 
late results. Similarly, one can develop several different theories of dynamic 
systems on time scales closely related to different numerical finite difference 
schemes. For instance, the standard delta calculus corresponds to forward 
(explicit) Euler scheme, the nable calculus corresponds to the implicit Euler 
scheme, and my proposition is related to the trapezoidal rule (and to the 
discrete gradient methods). Therefore there are no unique "natural" time 
scales analogues of dynamic systems. One can choose among many possi- 
bilities, including the above three approaches and the exact discretization 
(which explicitly exists only for a very limited class of differential equations). 
Another promising possibility is the so called "locally exact" discretization 
[13123]. 

The definitions presented in our paper seem to be entirely new as far as 
time scales are concerned but their discrete analogues (T = eZ) have been 
used since a long time. After completing this work I found a lot of references 
where rational or Cayley-like forms of the exponential function are used in 
the discrete case, see for instance |l6l[2ll[25l[26l|271|28l[29l[30]. It would 
be interesting to specify our results to the quantum calculus case [HI] , where 
the approach presented in this paper probably has not been applied yet. 

It is obvious that the proposed modification of the basic definitions should 
have an essential influence on many branches of the time scales calculus, 
including the theory of dynamic equations [U |5] , Hamiltonian systems [32] , 
and the Fourier and Laplace transforms [1]. 

Finally, we notice that exponential functions on time scales are defined 
(here and in other papers) on real time scales. It would be important to 
extend these definitions on the complex domain which is so natural for con- 
tinuous exponential functions. Such extensions are well known in the discrete 
case and in the quantum calculus, see for instance [211 [251 |29l 130] . 
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